Real-time dynamics of particle-hole excitations in Mott insulator-metal junctions. 



Luis G. G. V. Dias da Silva, 1,2 Khaled A. Al-Hassanieh, 3 Adrian 
E. Feiguin, Fernando A. Reboredo, 2 and Elbio Dagotto 1 ' 2 

1 Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996 
2 Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831 
3 Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 8754-5 
^Department of Physics and Astronomy, University of Wyoming, Laramie, Wyoming 82071, USA 

(Dated: March 23, 2010) 

Charge excitations in Mott insulators (Mis) are distinct Irom their band-insulator counterparts 
and they can provide a mechanism for energy harvesting in solar cells based on strongly correlated 
electronic materials. In this paper, we study the real-time dynamics of holon-doublon pairs in a MI 
connected to metallic leads using the time-dependent density matrix renormalization group method. 
The transfer of charge across the Mi-metal interface is controlled by both the electron-electron 
interaction strength within the MI and the voltage difference between the leads. We find an overall 
enhancement of the charge transfer as compared to the case of a (noninteracting) band insulator- 
metal interface with a similar band gap. Moreover, the propagation of holon-doublon excitations 
within the MI dynamically changes the spin-spin correlations, introducing time-dependent phase 
shifts in the spin structure factor. 

PACS numbers: 71.10.Fd, 71.35.-y 



I. INTRODUCTION 

The dynamics of excitations in strongly correlated elec- 
tronic materials (SCEMs) (Ref. U]) has been the subject 
of intense study in recent years. In a vast class of SCEMs, 
the ground state is a Mott insulator (MI), characterized 
by strong on-site repulsive interactions and a charge gap 
in the density of states. A paradigmatic model describ- 
ing this behavior is the one-dimensional (ID) Hubbard 
model at half-filling, in which elementary charge excita- 
tions involve either doubly occupied (doublons) or empty 
(holons) electronic sites. A doublon-holon pair is thus a 
charge neutral excitation, which, depending on the rela- 
tive strength of the on-site and nearest-site repulsion, can 
either form a bound state (a "Hubbard exciton") or re- 
main decoupled. These excitonic states in Hubbard-type 
systems have attracted considerable attention involving 
both theoretical^] and experimental^ investigations. 

A particular framework where the behavior of doubly 
occupied site excitations has acquired considerable rele- 
vance is in experiments involving cold atomic fermionic 
mixtures in optical latticesP While these experiments 
are performed under controlled conditions and allow for 
a substantial degree of tunability, the ability to prop- 
erly probe and characterize these systems is limited. A 
recently proposed technique to probe the spectrum of 
a Mott insulating state consists of dynamically creat- 
ing double occupancies by modulating the lattice depth 
with a time-dependent optical potential.^ When the lat- 
tice depth exceeds the Mott gap, double occupancy be- 
comes favorable and, by optical methods, it is po ssible 
to determine the appearance of the gapped modePHUl 
In addition, lattice-modulation spectroscopy has been 
proposed as a method to detect antifcrromagnetic or- 
dering and probe the nature (coherent or incoherent) of 
quasiparticle excitations.^ The decay mechanisms asso- 



ciated to the doublon excitations bring up fundamen- 
tal questions about thermalization and non-equilibrium 
dynamics EMU] 

In addition, technological applications may potentially 
arise from exploiting charge excitations in Mott insula- 
tors as a way to devise SCEM-based solar cells, which 
can offer structural and optical advantages over cur- 
rent devices made with semiconductor materials. One- 
dimensional Mott insulators (such as Sr 2 Cu03) are par- 
ticularly known for strong nonlinear optical response 
effectfPE3 that can be exploited in this context. The 
efficiency of SCEM-based solar cells will depend on sev- 
eral factors, including the performance of Mott insulator- 
metal junctions where the photocurrent will be gener- 
ated. A crucial question is whether charge excitations 
in the MI will be able to properly transfer the charge 
into the metallic contacts, thus establishing a steady- 
state photocurrent. 

Several questions arise in these regards: (i) Are the 
charge excitations long-lived? In other words, what are 
the effective decaying channels for the holon-doublon ex- 
citations into spin excitations inside the MI region? (ii) 
What are the effects of the dynamics of these excita- 
tions on the original correlations in the MI region? (iii) 
Can these charge excitations propagate across an inter- 
face with a non-correlated material? 

Point (i) was addressed by some of us in Ref. HH where 
the real-time dynamics of a holon-doublon pair was stud- 
ied in a ID Mott insulator. In that previous effort it 
was found that the decay to spin excitations in the un- 
derlying spin background is inefficient, and the pair is 
long-lived. The weak decay to spin-only excitations is 
particularly telling since the low-lying excitations for the 
ID Hubbard model at half-filling carry spin (spinon) and 
charge (holons) quantum numbers separately and propa- 
gate with different characteristic velocities, as described 
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FIG. 1: Schematic representation of the Ml+leads system 
studied here. At time r = 0, a doublon-holon pair is created 
at the center of the MI region. This pair propagates and 
eventually reaches the boundary interfaces. 



by the Tomonaga-Luttinger liquid picture. The ensu- 
ing spin-charge separation of these low-energy modes has 
been long studied and t he rea l-time dynamics has been 
explored in recent efforts (ASHSl 

Points (ii) and (iii) will be addressed in this paper. 
Here, we study the charge and spin dynamics across a 
Ml/metal junction, as doublon-holon excitations are cre- 
ated within the MI. A possible mechanism for such ex- 
citations arises from optical absorption at energies on 
the order of the Mott gap, producing exciton-like states. 
We study the time evolution of a localized doublon-holon 
pair initially created within the MI region. We consider 
a regime with weak nearest-neighbor electron-electron 
repulsion relative to the kinetic energy within the MI 
(or, equivalently, weak holon-doublon attraction), lead- 
ing to a spatial dissociation of the doublon-holon pair, 
with both excitations eventually reaching the Mi-metal 
boundaries. 

The main results are summarized as follows. A strong 
confinement of the holon-doublon pairs within the Ml 
region occurs for moderate to high values of the on-site 
interaction U. Nevertheless, charge transmission through 
the Mi-metal boundary can be favored by properly ad- 
justing the voltage difference between the metallic leads. 
These dynamical effects have a strong influence on the 
spin-spin correlations within the Ml region, effectively re- 
ducing the power-law decaying antiferromagnetic (AFM) 
order and introducing phase shifts in the spin structure 
factor. 

The paper is organized as follows. The model and the 
details of the calculation are described in Sec. HU The 
main results are presented in Sec. |III[ where we discuss 
the charge, double occupation, and spin-spin correlation 
dynamics (Sec. Ill A[ ), as well as compare results against 
noninteracting cases (Sec. IIIBl. The phase sh ifts in the 
spin structure factor are also discussed (Sec. |III C|). A 



II. MODEL AND METHODS 

The metal-MI-metal junction studied here is described 
by a ID Hubbard model with on-site and nearest- 
neighbor Coulomb interaction terms (representing the 
MI region) connected to noninteracting sites (represent- 
ing the metallic leads). The Hamiltonian reads H = 

HmI + Pleads + -^coupling, with: 



Hmi 

Pleads 
-^coupling 



a,i=N L + l ^ 

+v (m - i) (n i+ i - i) + h.c. 

i€R,L a,ieR,L 

-f (4^+ h.c.) (i) 

<T,i=N L ;i=N L +N Ml 



summary of the main results is given in Sec. |IV| 



where the sum in -ff C oupiing has only two terms (i = N-^ 
and i = Nj, + Nmi)- As usual, c- CT (c i(T ) represents the 
creation (destruction) operator of an electron with spin 
projection a at site i. 

In the model above, Hmi describes the central interact- 
ing region, with U and V being, respectively, the on-site 
and nearest-neighbor Coulomb repulsion, t" is the hop- 
ping matrix element between sites in the MI region. The 
second term (-ffi ea ds) represents the noninteracting leads, 
where t is the tight-binding hopping amplitude (taken 
as the unit of energy hereafter) and \it is the on-site di- 
agonal energy at the sites forming the lead £. Finally, 
^coupling describes the coupling between the leads and 
the MI involving a hopping amplitude t' . 

The total length of the ID chain representing the en- 
tire system is Ni, + Nmi + ATr, where JVl, Nmi, and Nr 
correspond to the number of sites in the left lead, the 
Mott insulator, and the right lead, respectively (referred 
to as a "Nl-Nmi-N-r configuration"). The key parame- 
ters governing the dynamics are the on-site interaction 
U within the MI and the on-site energy difference be- 
tween the leads, given by A[x = (/il — A*r)/2, which 
can be controlled by applying different voltages to the 
leads (for isolated leads in equilibrium at a fixed density, 
A/z will be the chemical potential difference between the 
leads). In the following, we consider symmetric chains 
with N L = N R = 20 and N m = 10 (for a total of 50 sites 
for the system) and assume a strong lead-MI tunneling 
amplitude (£' = t), which minimizes reflections at the 
interface due to a hopping matrix mismatch (although 
reflections due to the fact that t" ^ t still play a role). 

To investigate the dynamics of the excitations and 
tunneling into the metallic leads, we will focus on the 
regime of weakly bound doublon-holon pairs with negli- 
gible recombination probability. Following Rcf. 18, these 
constraints can be met by choosing U and V such that 
U/t" > 8 and V/t" < 0.6. Thus, here we select the 
values t" = 0.5t, V = 0.3t, and U > 4 for our studies. 
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As expressed before, all parameters with units of energy 
hereafter are given in units of the hopping t. 

Figure [I] depicts the metal-MI- metal junction. The 
equilibrium ("r < 0") ground-state \ip)o is obtained from 
static density-matrix renormalization group (DMRG) 
(Refs. I21I22|) calculations. At time t = 0, a doublon- 
holon pair is suddenly created, by acting with holon (h\) 
and doublon (d\) creation operators on the state \ip)o' 

\i>(r = 0))=hldl +1 \iP)o (2) 

With h\ = (1/V2) J2cr Civi 1 - and 4 = 

(1/V2) E CT 4"i- 

The site where the holon is created is chosen as p = 
Al + Ami/2 so that the doublon- holon pair is initially 
localized at the center of the MI region. This choice, 
although not crucial, keeps the system symmetric under 
the application of a particle-hole transformation followed 
by a reflection through the middle bond. 

The real time dynamics is obtained by time-evolving 
the original state using the time-dependent DMRG 
methocPEl and obtaining \tp(r)) = e- l " T \tp(0)). We 
find that a Suzuki- Trotter decomposition with time steps 
St = 0.05 — 0.1 and keeping 200 states during the time 
evolution provides an adequate choice for time-evolved 
quantities up to times r ~ 40. Since the sudden creation 
of an exciton via Eq. [2] is restricted to only two sites in 
the original system, we find that, starting from a well 
converged initial state, increasing the number of states 
during the time evolution does not significantly alter the 
main results. 

To probe the effects of the doublon-holon dynamics on 
the properties of the system, we calculate the expecta- 
tion values of local operators Oi at site i at each time 
step. We define the change from the corresponding equi- 
librium value as SOi(r) = (Oi)(r) — (Oi) cq , where (...) 
indicates the expectation value using the time-evolved 
state \iP(t)) at time t while (Oi) cq is the value calcu- 
lated at equilibrium ("r < 0"). Typical site operators 
considered here are the on-site charge rii and double oc- 
cupation Di = d\di. 

We also define the "charge transfer" to the right metal- 
lic lead as 

An R (r, A/i) = ^WC 7 ") ~ ( n i)cq = ^2 Sn t( T ), ( 3 ) 

ien ieR 

where the sum runs over sites in the right lead only. This 
quantity keeps track of the time-integrated charge that 
is transported into the right lead as a result of the cre- 
ation of the holon-doublon excitation at r = 0. A small 
charge transfer would indicate strong confinement of the 
doublon-holon pair within the MI region. 

We should point out that, as defined in Eq. [2] the 
doublon is created to the "right" of the holon (i.e., a 
"left /right" holon-doublon pair), making the right lead 
the "doublon side" and thus justifying the choice of the 
right lead for the definition of the charge transfer in Eq. 




20 10 20 j 30 40 50 




FIG. 2: Charge (m) at site i versus time r for a chain with 
N = 50 sites (20-10-20 configuration), using U = 4 and (a) 
Afj. = and (b) A/Lt = -1. 

|3j This is only an arbitrary convention since the creation 
of "left/right" and "right/left" holon-doublon pairs by 
optical absorption should occur with equal probability. 
In fact, the effective total charge transfer to the right 
metallic lead should take into account the contribution 
from a "right /left" holon-doublon pair as well. Given the 
symmetries of the system described above, this contribu- 
tion will be given by — Atir(t, — A/z) calculated for the 
left /right holon-doublon pair. Note that, once this con- 
tribution is taken into account, the net current is nonzero 
only for A/i 7^ 

Additionally, we investigate the spin-spin correlations 
away from equilibrium. We calculate the dynamical spin 
structure factor S(q,r), defined as: 

5 ^r) = ^ J2 ^- h ^{S*Sl){r), (4) 
j,fceMi 

with j, k spanning the interacting region (J, k € mi). 
As usual, a peak in S(q) at q = ir signals quasi-long- 
range antiferromagnetic order, given in real space by 
(SjSj +k ) ~ (— l) k /r a with a being the critical exponent 
on the order parameter. 
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FIG. 3: Spin structure factor S(q, r) in the MI region for 
U = 4 and Afi = 0. 



III. RESULTS 



FIG. 4: Nonequilibrium charge difference Srii (a,b), double 
occupation 5Di on each site (c,d), and spin structure factor 
S(q, t) in the MI region (e,f) versus time. The results are 
for U — 4 with different values of Afi: Afi = (a,c,e) and 
Afi = — 1 (b,d,f). Note the pronounced charge transfer for 
Afi = -I. 



A. Mi-leads charge transfer 

The typical dynamical behavior of the local charge 
(rii) (t) is depicted in Fig. [2} a for U = 4 and Afi = 0. 
At t < 0, the ground state of the system is particle-hole 
symmetric, with (rii) = 1 for all i. At r = 0, the doublon- 
holon pair is created, making (n p ) = 0, (n p+ i) — 2, and 
(D p ) = 0, (D p+1 ) = 1 (note that, while the total charge 
J2i n i 1S conserved, J2 i Di is not constant). Due to the 
on-site repulsion in the MI, these charge excitations move 
in opposite directions as a function of r, eventually reach- 
ing the Mi-metal interfaces at a time scale r ~ t t (for 
the parameters in Fig.[2j r r ~ 5.5). 

For Afj, ^ 0, the charge is expected to be non-uniformly 
distributed, with clear charge "plateaus" in each of the 
three regions. This is clearly depicted in Fig. [2}b, with 
MR > Ml (Afi < 0). In this case, the equilibrium charge 
distribution of the system is such that (n^i) > 1, and 
(njgij) < 1 while (n ie Mi) ~ 1 in the MI region. As the 
doublon excitation approaches the Mi-leads boundary, it 
is then partially transmitted to the right lead while, by 
symmetry, the holon excitation is partially transmitted 
to the left lead. 

An intuitive picture for such a charge transmis- 
sion/reflection can be formulated in terms of the single- 
particle density of states. The energy cost for the for- 
mation of the doublon is on the order of the Mott gap, 
Ami- Therefore, from purely energetic considerations, 
one expects an enhanced charge transport if there are 
enough states available in the leads at energies on the 



order of ± Ami/2 relatively to the Fermi energy. Static 
DMRG calculations for the equilibrium density of states 
(not shown) give Mott gap values Ami ~ 2 for U = 4 
and Ami ~ 8 for U = 10, which is on the order or larger 
than the typical half bandwidth ~ 2t of the noninter- 
acting chain. Therefore, for a fixed Afj,, one expects a 
decreasing charge transfer with increasing U as there are 
fewer states in the leads available at energies ~ ±Ami- 
On the other hand, a nonzero Afi can increase the charge 
transfer by allowing a "matching" of states available for 
\Afi\ - Ami/2. 

In fact, the reflection at the boundary can be reduced 
by changing Afi, as shown in Fig. [2]-b. In this case, the 
Fermi energy in the leads is aligned with the band, cre- 
ating states available for the excitation to "leak" into the 
leads. 

This simplified picture, however, does not take into 
account the existing correlation effects within the Mott 
insulator region. Additional insight on the effect of the 
doublon-holon dynamics can be provided by the spin-spin 
correlations within the MI region, via the spin structure 
factor. Fig. [3] shows the time evolution for S(q,r) cal- 
culated within the interacting (MI) region for U = 4. 
In equilibrium ( "r < 0" , shown as the first curve) , a 
pronounced peak is observed at q — tt in both cases, 
signaling AFM correlations. 

At intermediate times, this peak acquires "shoulders" 
which become more prominent around r ~ T r /2, which 
coincides with the time scale for which the doublon-holon 
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FIG. 5: Charge difference <5n;(r) on sites near the Mi-metal 
boundary as a function of time for U — 4 and (a) A/t = and 
(b)A/i = -l. 



pair is spatially separated but has not yet reached the 
boundary. For longer times, the peak is broadened, indi- 
cating a weaker AFM order within the MI region. The 
shoulders in S(q, r) are consistent with the formation of 
anti-phase domain walls in the MI region due to the spa- 
tial separation of the doublon and the holon, as we dis- 
cuss further in Sec. MI CI The fact that these features 
in the spin structure factor closely track the charge dy- 
namics indicates that charge excitations (rather than col- 
lective spin excitations) are the dominant decay chan- 
nel for the doublon-holon pair, in accordance to previous 
findings™ 

Figure [4] presents a visual summary of the dynamics 
of the charge, double occupation, and spin-spin correla- 
tions for U — 4 and different values of A/i. Figs. [4]a-b 
show (5nj(r) at all sites for A/i = and A/i = —1, respec- 
tively. An enhanced charge transmission into the leads 
is clearly seen for A/i = —1, consistent with the argu- 
ments given above. This additional charge is an effect of 
the doublon excitation tunneling into the leads. This is 
confirmed by analyzing the change in double occupation 
from equilibrium dD^r) (Fig. |4]c-d). Note that the dou- 
blon is "long lived" , as evidenced by the dynamics of the 
double occupation within the MI region as it reflects off 
the boundary. This is a consequence of the weak decay 
of the doublon for large U/t", as discussed in Ref. [T5] 

Additionally, there are interesting effects in the spin- 
spin correlation, as shown in the S(q,r) plots in Figs. [4] 
e-f. A clear splitting of the peak occurs as the doublon 
and holon excitations become separated within the MI, 
also seen in Fig. [3] More interestingly is the fact that this 
feature is enhanced for A/i < (Fig. |I]-f), a consequence 
of having a stronger charge tunneling between the MI 
and the l eads (i.e. stronger "doping" , as discussed in 

Sec. Imcl. 




A closer look on the charge dynamics at the interface 
is presented in Fig. [5j which depicts the behavior of the 
added charge Sn^r) at sites on both sides of the inter- 



FIG. 6: Charge transfer (Eq. [3} for (a) U = 4 and different 
values of A/i and (b) A/i = and different values of U, with 
t" — 0.5 in all cases. 



face. A weak odd-even modulation is present so only 
odd-numbered sites are shown for a meaningful compar- 
ison. The amplitude of the charge excitation decays as 
the "charge front" approaches the boundary. In addition, 
a reflection at the interface takes place for r ~ r r = 5.5 
in both cases. 

We can estimate the Mi-metal charge reflection by 
comparing the maximum amplitudes of Sn^r) at sites 
on each side of the interface. For A/i = (Fig. I5J- 
a), the maximum in Srii drops from 0.315 for i = 29 
(within the MI) to 0.071 for (in the leads), giving 

ife(max Sn 2 9 — max <5n 31 ) /max (5^29=0.7737. The rela- 
tive drop is smaller for A/i= —1 (Fig.pH-b), with i?«0.51. 

For longer times, further scattering of the doublon and 
holon excitations off the Mi-metal interface result in ad- 
ditional charge tunneling into the leads. The integrated 
charge transferred to the right lead (Atir, defined in Eq. 
[3| is plotted in Fig. [6}a for different values of A/i. For 
longer times, A«R(r) reaches an approximate plateau, 
indicating a steady state of the charge in the lead (for 
A/i = — 1.5 the plateau has not been fully reached on the 
maximum time used in the calculations). As expected 
from our previous results, the height of the plateau is 
larger for negative values of A/i. For a stronger inter- 
action U (larger Ami), the charge transfer is strongly 
suppressed, as shown in Fig. [6}b. 

These results indicate that the stronger the interaction 
U in the MI region, the more confined the holon-doublon 
pair becomes. This is clearly seen in Figs. [7ja,b) and 
(c,d), which depicts contour plots of Sn^r) and 8Di{r) 
for U — 4 and U = 10, respectively. For U = 10, most 
of the charge remains confined in the MI region, reflect- 
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FIG. 7: Nonequilibrium charge 5nj(r) (a,b) and double oc- 
cupation SDi(r) (c,d) on each site as a function of time for 
(7 = 4 (a,c) and U — 10 (b,d). Panels (e,f) show the corre- 
sponding spin structure factor S(q, r) in the MI region. 



ing off the interface in a periodic pattern. Interestingly, 
the long-lived holon and doublon excitations also "repel" 
each other, as seen for r ~ 2r r . 

Not surprisingly, increasing the Mott gap has a pro- 
nounced effect on the spin structure factor S(q, r). Figs. 
[7}e and f?\{ show S(q, r) for U = 4 and U = 10, respec- 
tively. For U = 4, the three-peak structure ("shoulders") 
disappear for r > r r , where a broad peak takes over. This 
is about the time scale for which the charge/spin fronts 
reach the boundary (Fig. [7]-b). For U — 10, the shoul- 
ders reappear at later times, consistent with the "beat- 
ing" in the doublon/holon pair as they reflect back and 
forth within the MI region and with each other. This is 
coupled with a much stronger confinement of the holon- 
doublon pair within the MI region, as seen in Figs. [7}b. 



B. Comparison with the non-interacting case 

An important question is how the results for the MI- 
leads charge transfer compare to a case where no electron- 
electron interactions are present in the central region. To 
address this question, we present results for two distinct 
cases: (i) a full "metallic" system, i.e., U = V = and 
t" = 0.5 in the Hmi term in Eq. [I] and (ii) a "band 
insulator" region, in which Hmi is replaced by 

N L +N MI 

H BI = -t" (l + S(-lY)clc l+la + h.c. (5) 

a,i=N L + l 




20 



10 20 : 30 40 50 



FIG. 8: Charge (rii) on site i versus time r for N = 50 sites, 
a non-interacting central region((7 = 0), and Afi — 0. (a) is 
using a Hamiltonian as in Eq. m with t" = 0.5 while in (b) 
Hbi replaces -Hmi an d t" = 1, 5 = 0.6 (band insulator). 



in the Hamiltonian. Notice that Hbi represents a Peierls 
chain, which has a charge gap for <5 ^ and the gap 
size will depend on t" and S. In the following, we choose 
t" = 1 and 5 = 0.6 so that the band gap is nearly the 
same as the Mott gap for the case shown in Fig. [2]-a 
(U = 4, V = 0.3, and *" = 0.5). 

Figures [8]-a and b show the occupation (n*) on the 
chain sites for the fully metallic and band insulator cases, 
respectively. The first case corresponds to the nonintcr- 
acting limit of the results presented in the previous sec- 
tion, with a fully metallic chain. There are still strong 
reflections at the boundary due to the hopping mismatch 
it" ^ t). More importantly, the holon-doublon charge 
excitations produce Friedel-type oscillations along the 
chain, forming a clear charge interference pattern over 
time (see contour plot in Fig.[8]-a). 

For the case where the central region is a band in- 
sulator, charge oscillations also occur, leading to a 
"checkerboard" pattern of alternating positive and neg- 
ative charges in the central sites (Fig. |8]-b). In this case, 
propagation of a "charge wavefront" is clearly suppressed 
as compared to the fully metallic case. 
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FIG. 9: Charge transfer comparison between the cases of 
interacting ((7=4) and noninteracting (!7=0) central regions. 
Results are also shown for a band insulator (Peierls chain) 
with a gap very similar to that of the Mott insulator with 
(7=4. (a,b) Ariij(r) versus r for A/i = and A/i = —1, 
respectively, (c) Time-averaged (steady-state) charge transfer 
(Anfl) T versus A/i for the MI and BI cases. 



The charge transfer to the right lead in both situations 
is clearly distinct from the interacting case. Figs. [9]-a,b 
show a comparison of the charge transfer at U = 4 (same 
as depicted in Fig. [6]) with results at U = for A/i = 
and A/z = — 1, respectively. 

The charge transfer is very limited for A/i = in 
the noninteracting cases, having essentially a zero time- 
average for the band insulator. By contrast, the case 
U = 4 shows a positive transfer even for A/i = due 
to the on-site repulsion within the central region. Note 
that the charge gap in the central region is nearly the 
same for both the Mott and band insulator cases, while 
the reflection at the boundary is much more accentuated 
for the latter. This highlights the role of the electron- 
electron interactions in the additional charge transfer for 
the U = 4 case. For A/i = — 1 (Fig. |9]4>) the charge 
transfer is improved in the noninteracting cases, as ex- 
pected. Overall, the charge transfer is still more effective 
in the case where the central region is a Mott insulator, 
as compared to the noninteracting cases. 

In fact, this holds for a wide range of values of A/z, 
as it can be seen in Fig. [9}c, where a time-averaged 
charge transfer (Atir) t (calculated by taking averages 
of Anjj(r) over a time interval Ar ~ 10 at longer times, 
away from the initial transient) is plotted as a function 
of A/i. In the MI case, (Anjj) T shows a broad peak, over 
a wide range of A/i (larger than the typical metallic half- 
bandwidth, ~ 2 in the units used), showing a positive 
charge transfer to the right lead even for positive values 
of A/i (i.e., /ifl < Hl). This is in sharp contrast with the 
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FIG. 10: Equilibrium spin structure factor S(q) for a Ml-only 
chain with N sites (no leads) calculated for states with charge 
N — Nh (Nh is the number of holes) and spin S z . Parameters 
are U = 4 and t" = 0.1. 



band insulator case, for which (Anu) T < for fin < fJ-L- 
Notice that, for A/i = 0, (Anjj) T > in the MI case 
while (An^) T w in the BI case (Fig. [j|c). The latter 
is an expected result: in a noninteracting system with 
an initially uniform charge distribution, an electron-hole 
excitation will produce charge oscillations about the ini- 
tial charge value (see Fig. |8| but, on average, no charge 
transfer takes place. 



C. Phase shifts in S(q) 

To understand the three-peak feature in the spin struc- 
ture factor (e.g., Fig. [3]), we have calculated S(q) for a 
Hubbard chain with N-^ = = and Nmi = 40 sites 
using static DMRG. We performed several DMRG cal- 
culations targeting states with different number of elec- 
trons N and total spin projection S z . More specifically, 
we considered half-filling (N — Nmi) and hole-doped 
(JV = /Vmi — Nh where iVh is the number of holes) cases, 
with S z = 0, 1/2, and calculated the corresponding S(q) 
for the different situations. 

Results are presented in Fig. 
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S z ~0, the state corresponds to the antiferromagnetically 
ordered ground-state of the system at half-filling, with 
a peak at q=ir. We then targeted states with different 
values of JVh and S z . This corresponds to calculating the 
ground-state of a system in which holes and/or spin flips 
are added. For iVj, = l, S z = l/2 (one hole) and N^ = 2, 
S z =0 (two holes) a double peak structure appears. 

This is reasonable since holons will introduce anti- 
phase domain walls (ADW) (i.e. "7r-shifts") in the AFM 
order and peaks at n ± nv^ are expected (where is 
the hole density). As spin-flips are added (for instance, 
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iVj, = 1, S z = 3/2) the double-peak becomes four peaks 
(not shown). In light of these static calculations, the 
appearance of "shoulders" in S(q, t) for r < r r (where 
r r is the time it takes for the excitations to reach the 
boundaries) can be accounted for as follows. 

The creation of the holon/doublon pairs at r = ef- 
fectively removes two magnetic moments from the MI 
region, making S*, +1 ^ = at the holon (doublon) site. 
This explains the decrease in the area under the S(q,r) 

curve at r = as / S(q,r)dq cx J2ieui{( S i ) 2 )( T )- For 
small r, S(q, r) still retains a single-peak structure ini- 
tially since holon and doublon are on adjacent sites and 
the corresponding phase shifts from ADWs cancel out. 

For r < r r , the MI region can be divided into two 
parts with magnetically distinct characteristics: (i) an 
"undoped" region where neither the doublon or holon ex- 
citations have arrived and which still retains AFM order 
and (ii) a "doped" region, which has been already "cov- 
ered" by either the doublon or the holon excitations. Part 

(i) contributes to a peak at q — 7r in S(q, r) while part 

(ii) contributes to the two shoulders, as we expect from 
the static calculations shown in Fig. |10| with = 2. 

For r ~ T r , the excitation (and the corresponding 
ADW) reaches the boundaries, making the phase shifts 
to cancel out again, suppressing the side peaks. As the 
holon and doublon reflect off the boundaries, the side 
peaks reappear up to time scales of order r ~ 2r r . At 
this time scale, doublon and holon excitations are again 
at nearest-neighbor sites and the respective ADW phase- 
shifts cancel once again, leading to a single-peak struc- 
ture. The process then repeats up to time scales on the 
order of the doublon decay time. Thus, this cycle is ob- 
served more clearly for larger values of U (e.g., Fig. [7|f) 
for which the doublon decay time is large enough. No- 
tice that S(q,T ~ 2r r ), although featuring a single peak 
at q = 7r, has a much smaller area than S(q, r = 0), in- 
dicating that the propagation of the holon-doublon pair 
significantly modifies the AFM correlations. 

IV. SUMMARY 

In summary, we have studied the real-time propagation 
of doublon-holon excitations in a Mott insulator (MI) 
connected to metallic leads. We analyze the dynamics 
of charge, double occupation, and spin-spin correlations 
within the MI, as well as the Mi-leads charge transfer. 

Our results indicate that the sharp change in the 
Hamiltonian at the Mi-metal interface hinders the charge 
transfer, suggesting that metals that closely resemble 
the structure of the Ml (Ref. [16]) would be required as 
charge collectors instead of standard doped semiconduc- 
tors. More specifically, we find that the charge transfer 
across the Mi-metal boundary is quite sensitive to micro- 
scopic parameters in the MI region, particularly the on- 
site interaction U. While U needs to be sufficiently large 
so that the doublon decay-time is larger than the typical 
time scale it takes to reach the boundary, doublon-holon 



tunneling into the leads is suppressed for very large val- 
ues of U. 

We believe two factors contribute to these results: (i) 
the increase in the Mott gap and (ii) the fact that, for 
large U, the nature of the doublon and holon excitations 
within the MI becomes more different than standard elec- 
tron and hole excitations in the noninteracting case. We 
have tested these hypothesis by comparing configurations 
with either a Mott insulator or a (noninteracting) band 
insulator in the central region: they show clear differ- 
ences in the charge transfer. 

The repulsive interaction within the Mott insulator 
region favors the transfer of the excess charge into the 
metallic leads even at A/j = 0. This is in sharp contrast 
with a noninteracting band insulator connected to leads 
with the same band-gap for which the net charge transfer 
is essentially zero at zero voltage difference. 

Noninteracting and interacting cases show clearly dis- 
tinct charge dynamics after the excitation. In the former 
case, Friedel-type charge oscillations in the central region 
are prominent at small times, while they are suppressed 
in the interacting case and a clear spatial separation be- 
tween hole-like and particle-like excitations occurs. 

Moreover, the propagation of holons and doublons 
within the MI region dynamically alters the AFM spin- 
spin correlations. In particular, extra "rr shifts" appear 
in the spin correlation functions as doublons and holons 
are spatially separated. We believe these qualitative find- 
ings will help on the prospect of making future solar cell 
devices using strongly correlated materials. 

An interesting aspect that remains open is the effect 
of a finite temperature in our results, a difficult prob- 
lem considering that the numerical study of transport 
in correlated systems at finite temperatures remains a 
very challenging subject. The energy scales of spin and 
charge excitations in ID Mott insulators can be quite dif- 
ferent (~ t" for gapped charge excitations and ~ (t") 2 /U 
for gapless spin excitations) and finite temp eratures can 
affect each of these channels differently! 25 * 26 ! In our case, 
however, since the doublon-holon pairs decay very weakly 
into spinon modes and the dynamics is governed by 
charge excitations, we expect the results to hold as long 
as temperatures are small compared to the charge gap. 
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